Formation of pure two-electron triplet states in weakly coupled quantum dots 

attached to ferromagnetic leads 
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Weakly coupled quantum dots in the Pauli spin blockade regime are considered with respect to 
spin-dependent transport. By attaching one half-metallic and one non-magnetic lead, the Pauli spin 
blockade if formed by a pure triplet state with spin moment S z = 1 or —1. Furthermore, additional 
spin blockade regimes emerge because of full occupation in states with opposite spin to that of the 
half-metallic lead. 

PACS numbers: 72.25.-b, 73.23.-b, 73.63. Kv 
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Spin-dependent transport in mesoscopic system is now 
a well-established research fieldi which has attained lots 
of attraction recently, both experimentalljiSiSiliSiSii and 
theoretically. 8 : 9 : 10 : 11 One major reason is the possibilities 
of using, for example, quantum dots (QDs) in spin-qubit 
readout technologies^; however, also fundamental ques- 
tions concerning how spin coherence and interactions in- 
fluence the transport are of main interest. Recent studies 
on weakly coupled QDs 1314 ' 15 revealed blockaded trans- 
port due to a unit population of the two-electron triplet 
states, the so-called Pauli spin blockade (PSB), which 
has also been theoretically veri&edm&H 

In serially coupled DQDs, the PSB can be found when 
the energy separation of the single electron levels be- 
tween the two QDs exceeds intra-dot charging energy 
subtracted by the inter-dot charging energy. The thus 
imposed inversion asymmetry on the DQD, i.e. the 
level off-set between the single electron levels in the two 
QDs, creates a unit probability to populate the two- 
electron triplet states in one direction of the bias voltage, 
henceforth referred to as forward biases. Because of the 
Pauli exclusion principle, there cannot flow any electrons 
through the DQD once the triplet state is fully populates, 
hence, the systems enters the PSB. 

In this paper, this system tuned into the PSB regime, 
where it is investigated with respect to spin-dependent 
transport. Replacing one of the leads by a ferromagnet, 
does not lift the unit probability of the population in the 
triplet state for forward biases. For definiteness, let for- 
ward biases mean that the chemical potential of the left 
lead is higher than that of the right lead. Thus, elec- 
trons flow from the left to the right lead, via the double 
quantum dot (DQD). Further, assume that the left lead 
is ferromagnetic. The broken spin symmetry creates a 
difference of the population probabilities for the three 
components of the triplet state which generates a finite 
spin moment in the DQD in the same direction as the ma- 
jority spin in the ferromagnetic lead. Particularly, in the 
limiting case where the ferromagnetic lead is a half-metal 
with spin a [a =T,4), the triplet is found to be in the 
pure two-electron configuration \o~) a\o~)b with unit prob- 
ability, where the subscripts indicate QDa/b, where QD^ 



(QDs) is adjacent to the left (right) lead. This can be un- 
derstood since the triplet consists of three states, hence, 
there is always at least one of the three states available 
for occupation independently of the spin-polarisation in 
the leads. Since the sum of population probabilities for 
these states is unity in the blockade regime, the popula- 
tion is distributed among the states in the most probable 
combination. 

For reverse biases the system behaves somewhat differ- 
ently when one of the non-magnetic leads is replaced by a 
half-metal. In the non-magnetic caseji&il the current is 
mediated through transitions between the two lowest sin- 
glet states and the lowest (spin-degenerate) one-electron 
state. The probabilities for these transitions are about 
1/2 to both leads. In the magnetic case with the left lead 
half- metallic (right lead non-magnetic), however, only 
electrons of the same spin as those in the half-metallic 
left lead are allowed to proceed through the DQD. Since 
both spin directions are available in the right lead and 
electrons of both spins are equally probable to tunnel 
into the DQD, the flow of electrons from the right tends 
to accumulate electrons in the state \ct)a\<j)b, that is, 
opposite to the spin in the half-metal. This accumula- 
tion is provided by tunneling processes like, for example, 
| initial) = \0)a\&)b -> |ct)a|0)b -> \o)a\o)b = | final), 
where it is understood that the electrons flow from the 
right to the left lead. 

If, instead, the right lead is half-metallic (left lead 
non-magnetic), this configuration also causes a block- 
aded current for forward biases. However, in this case 
the blockade is caused by a unit population probability of 
the triplet state \u)a\^)b, that is, the one with opposite 
spin to that of the half-metal. The reason is clear since 
electrons with the same spin as the half-metal can tun- 
nel through the DQD whereas those of opposite spin are 
blockaded by the absence of such spin states in the right 
lead. For reverse biases, the flowing electrons only have 
one spin direction (that of the half- metallic right lead). 
Since the singlet states require two electrons with oppo- 
site spin, these are unavailable for population, hence, for 
transport. In addition, since the probability for popula- 
tion in the triplet state is negligible when the left lead 
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is non-magnetic, 17 - only the one-electron state with the 
same spin as the half-metallic right lead is populated. 
This too leads to a suppressed current. 

The discussion is illustrated by considering two single 
level QDs (ea,£b, spin-degenerate) with intradot charg- 
ing energies (Ua = Ub = U), which are coupled by inter- 
dot charging ({/') and tunnelling (t) interactions, or more 

specincalh /WV H DQD = E<=A,n(£«r e i4<k + 
Uin^rni) + U'{n A ] + n A i){n B ] + n B ;) + J2 a ( td Aa d B<y + 
H.c). Here d Aa ^ Ba (dAa/Ba) creates (annihilates) an 
electron in QD a /B at the energy Eag/Bo^ while nAa = 
d Aa d,A<j and analogous for ubo- The PSB regime is ob- 
tained for ea — £b — Ae > 0, Ae < 2(/i — ea) < 5Ae, 
where /i is the equilibrium chemical potential, U = 2U' = 
2Ae, and weakly coupled QDs, e.g. £ = 2t/As <c 1. 
The condition Ae > imposes the required inversion 
asymmetry on the system (Ae < provides the PSB 
for reversed biases). As shown in Ref. 0, the range of 
chemical potentials satisfying 2(/i — sa) S (Ae, 5Ae) pro- 
vides a regime where the transport is mediated between 



one- and two-electron states in which the PSB is avail- 
able. Generally also low temperatures are required, e.g. 
ksT/U < 0.025. For a more detailed discussion on the 
mechanism behind the PSB regime for this system and 
the choice of parameters, as well as for the model of the 
DQD coupled to leads, the reader is referred to Ref. 

This configuration of the DQD and biases voltages V € 
(— 1, l)U/e, enables the analysis performed in Ref. ^3 
(to which the reader is referred for the general equations 
and notation) . Briefly a density matrix approach for the 
population number probabilities of the many-body eigen- 
states have been used. The DQD has 16 eigenstates with 
corresponding population probabilities P/vn, denoting 
the nth population probability of the TV-electron config- 
uration [N = 0, . . . , 4). In the Markovian approximation 
(which is sufficient for stationary processes) the equations 
for the population probabilities Pjy n , to the first order in 
the couplings T% R = 2irJ2 k eL/R \vka\ 2 5{uJ - E klJ ) to the 
left/right (L/R) lead, can be written as 



dl PN 



\—L,R ^ m 

- / jr^ w> ^ +lm [/^(Ajy- + i mi jv n )Pjv- n - /~(Ajv+l m ,iVn)-PjV+lm] J = 0, 
m 

N = 0,...,4, 



(1) 



where P_i„ = P 5n = 0. Here, A N+lm ^ Nn de- 
note the energies for the transitions \N,n)(N + l,m|, 



while T 



L/R 

Nn,N+lm 



ia<t/b*)nn+i, where 
(dAcr/Ba)NN+i 1S ^ ne transition matrix element for an 
electron exiting QDa/b- Moreover, f x (v) = f(u> — fJ, x ) 
is the Fermi function at the chemical potential fi x of the 
lead x = L, R, whereas fZ{w) = 1 — /+ (lo). Effects from 
off-diagonal occupation numbers, which only appear in 
the second order (and higher) in the couplings, are ne- 
glected since these include off-diagonal transition matrix 
elements to the second order (and higher) which gener- 
ally are small for £ <C 1. 

Ferromagnetic leads are introduced through the 

R _ 



parametrised couplings^ T ^ ^ 



= r (l±Pi)/2 and rf x 



r (l ±Pr)/2. Here, Pl/r € [—1,1] accounts for the 
magnetisation in the left/right lead such that p x = 
corresponds to a non-magnetic lead while p x = 1 (—1) 
corresponds to a half-metallic lead with spin f (J,). 

In the given bias voltage regime, it turns out 
that the only non-vanishing population probabilities 17 
are P\ m , m — 1,2, for the one-electron states 
\l,m) = a\a m ) A \0)B + P\0) AWm) b , 01,2 =t,l, and 
P2n, n — l,---,5, for the two-electron states |2,n) = 
\a n ) A Wn) B , n = 1,2, a h2 =U, |2,3) = [| T)a| I 



)b + | 1)a\ T>b]/>/2, and \2,n) = A n [\ ]) A \ 1)b - \ I 
) A \ r) B ]/V2 + B n \ U) A \0) +C n \0) A \ U)b, n = 4,5. 
Here, the two-electron eigenenergies are related as E 2 4 < 
£<2n < £-25, n — 1,2,3. Here also, a, (3, A n , B n , and 
C n are coefficients, which depend on the internal param- 
eters (eA,£ B ,t,U' \U) of the DQD, for the eigenstates 
expanded in the Fock basis. 

First consider forward biases in the rangeA* eV € 
(0.1, \)U . Then, the non-equilibrium conditions yield 
/i(A 2n ,im) ~ 1 and fR(A 2n ,im) ~ 0, where f L /ii(x) = 
f(x — Hl/r) is the Fermi function at the chemical poten- 
tial Hl/r, whereas A 2nAm = E 2n - E lm is the energy 
for the transition |1, to)(2, n|. This leads to the following 
equations for the population probabilities: 



Pirn = 



E 5 

"^2 t^L I ( j \mn 1 2 



T,a,m=l T a\( d Aa)?2 n \ Plr 



E' = iT R \{dB„) 



mn 12 
12 I 



-, TO = 1,2, 



1,.. 



(2a) 
,5(2b) 



These equations are further simplified by using the tran- 
sition matrix elements \ (d Aa / b a) 12™ \ 2 m Table||| A some- 
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TABLE I: Matrix elements for the relevant transitions ex- 
pressed in terms of the eigenvector coefficients. Here, L 2 = 
{f3A n /V2 + aB n ) 2 and R\ = (aA n /V2 + f3C n ) 2 . 
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what tedious but straightforward calculation yields 
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Pn 

Pl2 

P22 
P23 

P2n 



1+PRf ^\P2U 



l+PL \(3 

1 - p L f 1 +pr\ 2 ( a N ~ 



1 - p R V 1 + PL ) \P 

2 



P21, 



1 ~ Pl 1 + Pr 

1 + PL 1 - PR 

1-PlI+Pr 



P21, 



I + PL I ~ PR 



P21, 



1 - p L 1 + p R ( a L 



1+plI-pr\P R 



P21 



(3a) 

(3b) 

(3c) 
(3d) 

(3e) 



where L n = (3A n /\/2 + aB n and R n = aA n /y/2 + pC n 
Conservation of charge here implies that 1 = ^1 Pim - 
J2i p 2n, which gives 



P21 



1 



l+PR 



1+PL 
2 

' 1 



I- PL , 1 + PR A - PL 



a 

I- PL 



I- PR 1 + PL V 1 - PR 
1- PL l+PR 
1 — PR 1 + PL 



2\ 1 



(4) 



Now, consider the case p\ = 1 and pr ■ — 0, which cor- 
responds to letting the left lead be half-metallic with spin 
I whereas the right is non-magnetic. In this case the pop- 
ulation number P 2 i reduces to P 2 i = + (a//3)*/2], 
where a 2 = £ 2 /[(l + ^/nK 2 ) 2 + £ 2 ] and /3 2 = (1 + 

C 2 )/[(l + v/l+^) 2 + e\ giving (a//?) 2 = £ 2 /(l + £ 2 )- 
Hence, 



P21 



i+e 

2 + 3£ 2 



1 as £ ^ 0. 



(5) 



Physically, this results mean that the population prob- 
ability of the triplet configuration | T)fl is unity. 
Consistency with Eq. is immediate, since all but Eq. 



(|3a|l vanish for p L = 1, while Pn oc (a//3) 2 -> 0, £ -> 0. 
The result is also consistent with the normal PSB, since 
the population probability of the triplet states must be 
unity. Hence, while only spin | electrons can tunnel from 
the left leads into the DQD. the only configuration avail- 
able for these electrons is | T)b- 

Similarly let pl = and pr = — 1, which corresponds 
to a non-magnetic left lead and half-metallic right lead 
with spin J,. Then, one finds that P21 = 1 apparently in- 
dependently of the coupling strength between the QDs. 
Thus, the state | t) ^1 1 T)b acquires a unit population 
probability for a spin J, half-metallic lead to the right. 
Consistency with Eq. © holds since all equations van- 
ish for Pr = — 1. A more careful analysis shows that 
the approximation, e.g. putting /L(A 2ri! i m ) — 1 and 
fR(A2n,im) = 0, provides a slightly overestimated popu- 
lation probability for | T)b- Nonetheless, configur- 
ing the system in this fashion indeed yields a significantly 
strong blockade. 

By turning to the reverse biasing case, e.g. eV £ 
(— 1, — 0A)U, one finds from analogous arguments, e.g. 
/i(A 2 „,i m ) « and /fl(A 2n ,i m ) « 1, that 
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PL 
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(7) 



ft follows from Eq. J7J) that a half- metallic spin | left lead 
and non-magnetic right lead, e.g. pl — I and pr = 0, 
yields P 22 = 1 independently of £ since all terms in the 
denominator but the first vanish for this configuration. 
Moreover, all relations in Eq. JBJ vanish for pl = 1. 
Physically this implies that the two-electron state | j 
)a\ 1)b acquires a unit population probability for low 
temperatures whenever the left lead is a spin | half- metal. 

Assuming instead that the right lead is a spin J, half- 
metal whereas the left lead is non-magnetic, e.g. pl = 
and pr = —1, leads to the conclusion that P 22 = 2£ 2 / (1 + 
3£ 2 ) — ► as f -> 0. Further, all relations in Eq. © but 
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FIG. 1: (Colour online) DQD system with a) pl = 1, pr = 0, 
b) pl ~ 0, pr = — 1. Population (solid) and spin-polarisation 
(dotted) of the DQD as function of the bias voltage, as well 
as the spin-polarisations of QDa (dash-dotted) and QDs 
(dashed). Here, £ = 0.01, k B T/U = 0.01. 

Eq. (|6b|) vanish in this case, whereas 

P 12 = -( P 22 = 1 + ^1 -» 1 as £ -> 0, (8) 
2Va/ l + 3^ 2 

that is, the one-electron state a\ |)a|0) + /^lO)^! 1)b ac- 
quires an almost full occupation for weakly coupled QDs. 
Hence, the DQD end up in a spin-polarized configuration, 
however, with only one instead of two electrons. 

The plots in Fig. ^ display examples of the population 
(solid) and spin-polarisation (dotted) in the DQD when 
coupled to one non-magnetic lead and one half-metallic 
lead a) p L = 1, p R = and b) p L = 0, p R = -1. 
Case a) shows regimes for both positive and negative bi- 
ases in which the DQD is populated by two electrons, 
where the spin-polarisation is 1 and — 1 respectively, cor- 
responding to the cases treated in the analytical calcu- 
lation above. Studying the individual spins of the QDs, 
QDj4 (dash-dotted) and QD# (dashed), shows that the 
two-electron regimes are indeed constituted by the triplet 
state | T)b an d | [)a\ 1)b i0r positive and nega- 

tive biases, respectively. In case b), there is only one 
two-electron regime in which the spin-polarisation is 1. 
Thus, the two-electron state consists of only the triplet 
state | ])a\ T)b f° r positive biases. For negative biases, 
the DQD is populated by a single spin | electron in the 
a \ I)a|0)s + /3|0)a| 1)b state, as expected from the an- 
alytical analysis. The weak coupling between the QDs, 
i.e. £ -C 1, implies that a as and j3 « 1, that is, the 
single spin J, electron is mainly located in QDg. 

Although the main focus in the plots shown in Fig. 
^is on the blockade regimes, there are a few details in 
that require some extra discussion. However, the struc- 
ture for biases |eV|/Ae > 2 goes beyond the focus of 
the present discussion, since this regime requires an ad- 
ditional investigation of the population numbers. In case 
a), the dip in the population around equilibrium (eV = 0) 
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FIG. 2: (Colour online) Magnitude of the current through 
the system for various strengths of the coupling £ of the QDs. 
Other parameters as in Fig. Q 



is due to that the lowest two-electron singlet states (as 
well as the triplet) are nearly degenerate with the lowest 
one-electron states, hence, the energy required for tran- 
sitions between those states is vanishingly small. This 
gives rise to an averaged population which is less than 
two electrons. The spin-polarization goes through zero at 
eV = since the lowest two-electron state is spin singlet, 
hence the spin-polarization vanishes. The high bias volt- 
age dips (eV/Ae « —1.6 and eV/Ae « 2) are due to the 
lifting of the triplet blockades caused by the availability 
of more transitions. In case b), the population goes from 
one to two electrons around equilibrium, which verifies 
the analytical discussion above. The spin-polarization of 
the DQD, however, has a regime for small positive bi- 
ases where the spin "f is slightly favored in both QDs but 
where the polarization has not saturated. This situation 
arises since the two electrons in the DQD is distributed 
among the three triplet states, although there is a slight 
overweight on the | ^)a\ T)b state. As is seen, this state 
is eventually fully occupied as the bias voltage becomes 
sufficiently large. 

The current through the system is calculated by the 
formula 22 J L = -(e//i)tr Im / T L [/l(w)G> (lj) + (1 - 

f L )(uj)G < (uj)}duj, where G</> = {G^ N+ln } Nmn is 
the matrix lesser /greater Green function (GF) of the 
DQD, whereas T L = £ CT T L a {\{d Aa )^ n N+1 \ 2 } Nmn is the 
matrix of the coupling to the left lead times the ma- 
trix elements for the transitions. The lesser and greater 
GFs are identified bjiii G^ m N+ln (u) = -i2TtP Nm S(oj - 
A N+ i ntNm ) and G^ mN+ln (uj) = i2nP N+ i n 5(uj - 
A N+ i„. Nm ), respectively. 

First consider case a), i.e. pi, = 1, p R = 0. Then, 
in the forward PSB regime, i.e. /i(A2 n ,i-m) = lj the 
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expression for the current reduces to 



Jl 



eT a 2 /2 
h l + (a//?) 2 /2 



0+, £ 



0. 



(9) 



Hence, for weakly coupled QDs the current is significantly 
suppressed, and decreases as £ 2 , £ — > 0. For negative 
biases the system enters the PSB regime as /r,(A 2ni i m ) = 
0, giving Jl = because of the simplifications leading to 
Eqs. JBJ| and J7J. Nonetheless, the result clearly indicates 
a significant current suppression. 

The simplified analysis is verified by the calculated cur- 
rent, displayed in Fig. EJa), showing the magnitude of the 
current as function of the bias voltage for various cou- 
pling strengths £. The plots clearly demonstrates that 
the forward current in the blockade regime decreases as 
£ 2 . In the reverse bias blockade, however, the current 
grows roughly by £ as £ — > 0. The reason is that elec- 
trons in the right lead enter the DQD via transitions 
between the triplet and one-electron states with proba- 
bility a 2 ~ £ 2 . Thus, for extremely weakly coupled QDs, 
the triplet state acquires a non unit population probabil- 
ity. The remaining probability is distributed among the 
singlet and one-electron states, which lifts the blockade. 
The pronounced y-shape of the current for reverse bi- 
ases arise from the fact that the population in the triplet 
state is not a completely flat function of the bias. In- 
stead, the state | J.)^| 1)b is maximally occupied around 
eV/Ae w — 1 which, hence, leads to the minimal current 
at this bias. The dip at equilibrium is because the log- 
arithm of the modulus of the current, i.e. log 10 \Jl\, is 
shown rather than the smooth function Jl- For higher 
bias voltages, e.g. |eV|/Ae > 2, the triplet blockade is 
lifted, which is shown as a significant increase in the cur- 
rent (several orders of magnitude) through the DQD. 

The case b), i.e. pl = 0, pn = — 1, is treated similarly. 
For positive the biases the above analysis provides Jl — 
0, which indicates that the current should be significantly 



suppressed in the forward blockade regime. For negative 
biases in the blockade regime, i.e. /z(A2n,im) = 0, the 
expression for the current reduces to 



er 



h l + 2(a//3) 2 



0-, £^0. 



(10) 



The plots in Fig. [21 b) display the calculated current, 
showing that the forward current is dramatically reduced 
in the blockade regime. The blockade is strengthened by 
a weakened QD coupling £, which is reasonable from the 
analysis in the non-magnetic case^i The reverse current, 
in the regime with a stable one-electron configuration, is 
also shown to be substantially suppressed (decreasing as 
£ 2 ), as expected from Eq. i|l(J|) . The V-shaped current 
for forward biases stems from the same reason as in case 
a). 

In the present analysis a few approximations have been 
introduced. First, the density matrix of population prob- 
abilities is diagonal. This is reasonable since the prob- 
ability for off-diagonal transitions, i.e. transitions be- 
tween different states with the same number of electrons 
vanish in the sequential tunnelling limit. Taking this 
limit is justified by the large off-set Ae > between the 
QD levelsiii The presented results are also confirmed in 
models with more levels in the QDs»22i Secondly, ferro- 
magnetic leads may induce a spin-split of the states in 

L /i? 

the DQD which is of the order of the couplings T a be- 
tween the DQD and the leads^ Weak couplings T„ , 
however, only marginally modifies the presented results. 

The conditions required in order to obtain the PSB 
regime are experimentally accessibleiiiiiiii 5 . Half-metallic 
leads may be constructed from NiMnSb^i which is com- 
patible with existing semiconductor technology. 25 An- 
other interesting half-metallic material is» Cr02 because 
of its favorable switching behaviour at low fieldsjSJ. 
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